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SUMMARY 


The properties of a Navier-Stokes solution of a shock-separated turbulent 
flow over a flat wall are investigated. Refinements of an algebraic relaxation 
turbulence model previously shown to be of value for the simulation of separated 
flows are presented. A simplified analysis applicable near an adiabatic wall 
is developed and used to help verify the accuracy of the numerical solution. 
Features of the time-dependent response of a turbulent boundary layer to shock 
impingement are presented. 


INTRODUCTION 


Computers now available are capable of practical calculations of complex 
flow fields, including separated turbulent boundary layers. However, as 
discussed in reference 1, development of adequate turbulence models is a pacing 
item that impedes progress toward that goal. Recent improvements in numerical 
methods, such as those described in reference 2, have made it feasible to test 
a variety of modifications of existing turbulence models (see, e.g., refer- 
ences 3-5 and the more comprehensive reference lists therein) . For engineering 
purposes, it would be practical to use simplified models calibrated from 
experiments conducted at nearby flow conditions. The shortage of experimental 
information on separated flows prevents calibration with precision at this time. 
Nevertheless, it seems worthwhile to proceed with the development of computer 
codes for complex flows based on simplified turbulence models that can be 
adjusted to accommodate the existing experiments. With presently available 
information, it may be possible to accomplish this for high Reynolds number 
flows at Mach numbers up to 3 using the boundary-layer approximation for 
viscous and Reynolds stress terms in the layers near solid surfaces. More 
complicated procedures can be incorporated when they are justified or can be 
used to aid in the calibration of the simplified models. 

In reference 5, it was shown that an inner layer algebraic eddy viscosity 
model used by Clauser (ref. 6) provides better agreement with two separated flow 
experiments than more conventional models. In this paper, the properties of a 
Navier-Stokes solution of a shock-separated flow based on a variation of that 
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model are Investigated. A simplified analysis applicable near an adiabatic 
wall is presented. Favorable comparisons of results from this analysis with 
the Navier-^Stokes' solution show that the botmdary-layer approximation is valid 
near the wall for this case and that the Navier-Stokes solution is accurate. 

As an example of the type of information that can result from such numerical 
solutions, features of the time-dependent response of a boundary layer to 
impingement of a shock wave are presented. 


SYMBOLS 


A constant (eq. (22)) 

B constant (eq. (23)) 

B_ constant in turbulence model (eq. (6)) 

r 3 

C constant (eq. (17)) 

Cp specific heat at constant pressure 

fj(y+) universal function (eq. (18)) 
universal function (eq. (19)) 
k Karman constant (0.4) 

K Clauser constant (0.016) 

p pressure, N/m2(lb/ft2) 

P turbulent Prandtl number (0.9) 

RT 

Sj, S^ constants in Sutherland viscosity law (eq. (20)) 
T temperature, K (®R) 

u velocity in x-direction, m/sec (ft/sec) 

u friction velocity (eq. (4)), m/sec (ft/sec) 

T2 

u friction velocity (eq. (5)), m/sec (ft/sec) 

t3 

V constant in Van Driest damping factor (eq. (8)) 

V„ Van Driest damping factor (eq. (8)) 

r 

x,y Cartesian coordinates, m (ft) 

y^ law of the wall variable (eq. (9)) 
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6 


boundary-layer thickness, m (ft) 

6* kinematic displacement thickness (eq. (10)), m (ft) 

X relaxation parameter (eq, (5)) 

y molecular coefficient of viscosity, kg/m-sec (slugs/ft-sec) 

turbulent eddy viscosity coefficient (eqs. (1),(2)) kg/m-sec 
(slugs/ft-sec) 

p gas density, kg/m^ (slugs/ft^) 

T shear stress (eq. (11)), N/m^ (Ib/ft^) 

Subscripts : 

0 initial profile at station ahead of interaction 

2 value of y at which u^^ is evaluated (eqs. (3), (4)) 

3 value of y at which evaluated (eq. (6)) 

1 position of inviscid shock impingement on wall 
max maximum velocity in profile 

w value at wall 


METHOD 


Numerical Method for Navier-Stokes Solutions 

The basic numerical method used in this investigation is described in 
reference 7. Recently, MacCormack has improved the method such that the calcu- 
lations require an order of magnitude less computation time than formerly (see 
ref. 2). The compressible Navier-Stokes equations to be solved are also 
listed in reference 7. 


Experiment Used for Comparison 

The experimental flow field (ref. 8) is depicted in figure 1. A shock wave 
generated by a plate set at 13® to the free stream Impinges on the boundary 
layer on the upper wind tunnel wall. A separation bubble containing reversed 
flow forms and extends upstream of the inviscid shock impingement point. A 
pattern containing induced and reflected shocks forms. The free-stream Mach 
number is 3 and the Reynolds number based on Initial boundary-layer thickness 
ahead of the interaction is 10^. 
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Turbulence Model 


The two-layer algebraic eddy viscosity model used in the present calcula- 
tions is defined by the following equations: 


f Pa ^ , (y 5 yj' 


("t). 




inner 


k = 0.4 


( 1 ) 


L. Pkyu^3V2 


• (y =■ ^2). 


/y \ = Kpu 6* , K = 0.0168 

V ^ max ’ 

outer 


( 2 ) 


where 1^2 evaluated at y = y, with determined from the 

relations 


' y=y, 


= 100 


“,2 • (■'y I?) 


y/ 

•' TT=' 


y=y. 


(3) 

(4) 


The quantity u^^ is determined from the relaxation formula 


du u - “ u „ 

t 3 _ T3eq ^ 

dx Ai5 Q 


X = 5 


(5) 


where 


^^T3eq 9y) 


y=y. 


, Bp3 = 1.18 


with yg determined from the relation 

3v 


^3 


' ^ y=y, 


= 2000 


(6) 


(7) 


The quantity V_ is the Van Driest damping factor 

r 


Vp = 1 - exp(-y+/V) , V = 18 


( 8 ) 


where 


+ '’2“t3y 

jT sr 


(9) 
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The transition from the inner to the outer formula takes place at the minimum 
value of y at which (U^-) inner ^ ^^t^outer* kinematic displacement thick- 

ness 6* is given by 


6 


* 



( 10 ) 


The foregoing model differs from that developed in reference 5 in several 
respects. Replacement of p with at y 72 equation (1) has 

the effect of removing the dependence of on variations of temperature 

and density in the viscous sublayer. Although future experimental data may 
show such a dependence, it seems preferable at this time to make dependent 

on alone, and thus to maintain a close correspondence with incompressible 

flows. Additionally, the factor equation (6) allows a close corre- 

spondence of the present inner layer Clauser model with more conventional models 
for flows with zero pressure gradient. Finally, use of y"^ as defined in 
equation (9) in the Van Driest damping factor instead of the definition used in 
reference 5 requires rescaling of V to V = 18. 


Simplified Analysis Applicable Near an Adiabatic Wall 

The following approximate relationships, obtained from the compressible 
Navier-Stokes equations, are useful for checking the adequacy of the mesh used 
for the niamerical solution. They may also be useful for deducing values of 
parameters in the turbulence model from experimental data obtained from sepa- 
rated and attached boundary-layer flows with pressure gradients. Upon neglecting 
the convection and inertia terms, the steady-state x-momentum and energy equa- 
tions can be approximated by 


X = (p + Pt) = -^2 + (y - y 2 > ^ 


C T C T 

= “p~ T “2 (adiabatic wall) 
RT- RT 


( 11 ) 

( 12 ) 


Replacement of the y-momentum equation with 9p/9y = 0 and use of the equation 
of state for a perfect gas yields 



With the definition 


du/ C 

I 1^2 

I -2- du/C 

V P2 


(y < y2> 
(y > y2> 


(13) 


(14) 
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and substitution of the foregoing inner-layer eddy viscosity model into 
tlon (11) , the following expression for u+ can be derived ; 

equa- 

where 

u+ = (1 - p+)f, (y+) + ^ f„(y+) 
i yj 2 

(15) 


+ _ dp_/dx _ , _ Iw 

^ ^2 T T 

^2 ^2 

(16) 


■ "t2 [■% * 

(17) 


y+ 

(18) 


4.. f y V 1 + 

1+V.yW) 

0 ^ 

(19) 

The approximations 
Vp - 1. 

for f j and f^ on the right apply at y"*" > 60 where 


To obtain the 
law in the form 

relation between u and u, we use the Sutherland viscosity 


SjT3/2 
~ $2 + T 

. (20) 

Linearization of the Sutherland relation and substitution of equations (12) 
and (13) into (14) leads to 




+ , 
u = i 

y 

^ (1 - i AB2 u 2) ^ (u < U 2 > 

+ .®V2, ra + Bu)(l-Bu2)] , ^ 

(21) 

where 

2 + aSj/T^ 

^ 1 + S 2 /T 

2 w 

(22) 


/2C T 2 
\ / 

(23) 
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RESULTS Al^D DISCUSSION 


A time-dependent solution of the compressible Navier-Stokes equations based 
on the foregoing eddy viscosity model has been carried out. The initial flow 
field is uniform in the x-direction with a boundary-layer profile corresponding 
to an upstream station in the experiment of reference 8. At the lower boundary 
of the computational field, boundary conditions are imposed corresponding to 
the shock wave in the experiment (see fig. 1). During the calculation, the shock 
wave grows toward the upper nozzle wall. Eventually, a steady state is reached 
corresponding to the conditions at which experimental measurements were made. 

As an illustration of the type of information that can result from such calcula- 
tions, features of the time-dependent response of the boundary layer to growth 
of the incident, induced, and reflected shock pattern will be presented. How- 
ever, the simplified analysis in the preceding section will first be used to 
check the adequacy of the computational mesh. 

The steady-state Navier-Stokes solution was used to compute dimensionless 
profiles of u***" versus y”^ according to equations (9) and (21) . Three such 
profiles are contained in figure 2. The circles represent the Navier-Stokes 
solution and are at the computational mesh points. The dashed lines are obtained 
from the simplified analysis (eq. 15)), using values of p"^ and y’t evaluated 
from the Navier-Stokes solution according to equations (3), (4), (9) and (16). 

The upper two profiles are at stations aft of reattachment where relatively 

small pressure gradients are present. The bottom profile is in the region of 

constant pressure ahead of the separated region, which extends from 

-2.6 < (x - x^/6q)'< 0.5. The close correspondence between the approximate 

and numerical solutions leads to two conclusions: (1) Use of the boundary- 

layer approximation with additional neglect of convection and inertia terms 

is a valid approximation near the wall for small pressure gradients; and 

(2) the computational mesh used for the Navier-Stokes solution provides adequate 

resolution. However, are the same conclusions valid at stations where strong 

pressure gradients exist? 

Profiles near separation and in the middle of the reversed flow bubble 
are shown in figure 3. Again the Navier-Stokes solution is represented by 
s 3 nnbols and the simplified boundary-layer approximation by dashed lines. 

Near separation, the approximate results are invalid for values of y"^ greater 
than about 300 because of neglect of the convection and inertia terms. In the 
middle of the separation (lower curve) , the simplified analysis retains validity 
to large values of y+. The resolution of the Navier-Stokes solution is again 
shown to be adequate. Profiles near reattachment are shown in figure 4. The 
two solutions agree closely over a large range of in this region. 

With confidence in the resolution of the Navier-Stokes solution established, 
it is of interest to observe the time-dependent response of the boundary layer 
to shock impingement. Figure 5 contains plots of wall pressure distribution 
at a series of time intervals after the start of the calculation. Shortly 
after the shock reaches the boundary layer and a reflected shock has formed, the 
wall pressure rise is steep and extends about one boundary-layer thickness 
upstream of the inviscid shock impingement point, which occurs at 0 on the 


1489 



abscissa scale. During succeeding time intervals, the pressure gradient 
decreases and the pressure rise moves upstream. Eventually, a steady state is 
reached in which the initial pressure rise occurs about three boundary-layer 
thicknesses ahead of the inviscid shock impingement point. The plot at the top 
includes the corresponding experimental steady-state pressure distribution from 
reference 8 for comparison. 

Figure 6 contains plots of the skin-friction distributions after the same 
series of time intervals. In the plot at the top, the calculations are in close 
agreement with a Preston tube measurement of the initial skin-friction coeffi- 
cient Cf and oil-flow observations of separation and reattachment points 
(ref. 8). Skin-friction measurements were not made at other stations in this 
experiment . 
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